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Abstract. We investigate a two-dimensional single-band Hubbard model with a nearest- 
neighbor hopping. We treat a commensurate collinear order as well as incommensurate spiral 
magnetic phases at a finite temperature using a Hubbard-Stratonovich transformation with 
O ' a two-field representation and solve this problem in a static approximation. We argue that 
• temperature dramatically influence the collinear and spiral magnetic phases, phase separation 

■ in the vicinity of half-filling. The results imply a possible interpretation of unusual behavior of 
! magnetic properties of single-layer cuprates. 

Introduction 

Zj i Investigation of two-dimensional (2D) electronic systems attracts substantial interest, which 

has been stimulated by the discovery of high-temperature superconducting cuprates. It is gen- 
CN ■ erally accepted that superconducting and magnetic propeties of cuprates are closely related, 

■ but magnetic properties are formed by Cw-layers and interlayer coupling is weak. At half- 
'. filling, the cuprates are antiferromagnetically ordered, evolution of their magnetic properties 

O 



with doping and temperature is an interesting challenge. For example, neutron scattering in 
^ _ La2-xSrxCu04 reveals coexistence of both commensurate and incommensurate magnetic struc- 
Q ■ tures in the vicinity of half-filling at low temperature, that pass to antiferromagnetic with rise 
of temperature [1, 2, 3, 4]. Recently the authors [5] considered the ground-state magnetic phase 
diagram of the two-dimensional single-band Hubbard model with nearest and next-nearest- 
neighbor hopping in terms of electronic density and interaction. They treated commensurate 
^ I ferromagnetic and antiferromagnetic as well as incommensurate spiral magnetic phases using 
■ the mean field (MF) approximation. First-order magnetic transitions with changing chemical 
potential, resulting in a phase separation (PS) in terms of density, was found between ferro- 
magnetic, antiferromagnetic, and spiral magnetic phases. 

Here phase diagram is investigated depending on temperature. Our calculations are based 
on a two-dimensional single-band Hubbard model. We use a Hubbard-Stratonovich transfor- 
mation with a two-field representation and solve this problem in a static approximation. It 
allows us, on the one hand, to obtain a solution that coincides with the MF approximation at 
zero temperature and, on the other hand, to investigate the temperature behavior of PS and 
phase transformation taking into account thermodynamic transverse spin fluctuations. 

First of all, we want to emphasize the peculiarity of the 2D magnetic systems. The series 
of works proved rather directly impossibility of existence of a spontaneous magnetic order and 
phase separation in a two-dimensional Hubbard model [6, 7, 8] (at finite temperature the results 
[6, 7] coincide with Mermin's theorem). At the same time, there are a lot of papers using the 2D 
Hubbard model for describing cuprates (see, for example, [9, 10, 11, 12] and references in them). 
Usually, it is believed that weak interlayer coupling stabilizes a magnetic order, which allows us 



to consider fluctuations around this broken-symmetry state even at finite temperature. There- 
fore, one can choose approximations which depress the 2D divergence of magnetic fluctuations 
(for example, a mean field approximation, a dynamic mean field approximation and so on). In 
addition, phase separation and inhomogeneous magnetic state in the 2D Hubbard model have 
been obtained recently by calculation methods using an extension of the dynamic mean-field 
theory and advances in a quantum Monte Carlo techniques [9, 10, 11]. We want to emphasize 
that exact results of [6, 7, 8] were received for homogeneous states; therefore, for inhomogeneous 
case the long wave fluctuations are cut off for the mean length of homogeneity. It reminds the 
stability of a 2D graphene layer where ripples appear due to thermal and quantum fluctuations 
and this stabilizes the crystal [!>>]. Though the problems of a 2D-crystal stability and of an 
inhomogeneous magnetic order in the 2D Hubbard model have not been solved yet, the results 
of the calculation methods give us an additional basis for using approximations which stabilize 
states with a broken symmetry. In this paper we neglect magnetic longitudinal fluctuations and 
thereby depress the 2D divergence. 



Method 

We consider the Hamiltonian of the Hubbard model on the square lattice H = Hq + Hint, 



E 



/ ^ ''t3'-'i,s'-'j,sy 
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(1) 



where j is a number of site, s is a spin projection, ti 



-t for the nearest-neighbor site 



and if not, c- 



is a creation (annihilation) electronic operator on the site j and U is 



the electronic (Hubbard) on-site interaction. We use a two-field formalism in the Hubbard- 
Stratonovich transformation, therefore more convenient representation of the on-site interaction 
in form 



H^nt = UY, = UJ2 N/4 - (S, ■ e 



(2) 



where is the (arbitrarily chosen) j-dependent unit vector and we introduce the site density 



Tin 



%.t + "^hi ^^"^ ^^^^ magnetization operator Sj = 1/2 s' cts^s.s'C, 



{<Ts,s' are the elements of vector cr = (a^, Oy, a^) formed by the Pauli spin matrices). 

Thermodynamical properties are defined by partition function of the grand canonical en- 
semble: 



Tr 



exp < —/3{H — yun 



Tr 



exp |-/3(-ffo - M)| TV exp < - / Hint 



'j)dT 



■ (3) 



Here /x is a chemical potential, is a time ordering operator, Tr is summation over a complete 
set of states, Hintir) = exp (— r_ffo)-f^mt exp (rf^o) is an operator in the interaction represen- 
tation, /3 = l/Zc^T, T is temperature. Partition function can be rewrite with the help of a 
Hubbard-Stratonovich transformation 



Z = I n^v,(r)DO(r)^[v,C], 



(4) 



Z[v,C]=exp{-/3F(v,C)} = exp|-f/^ J rfr[v2(r) + C|(r)] U*[v, C], 



Z*[^r,C] = Tr 



exp \-l3{Ho 



p 



V;,-(v, C, r) = f/ zO(r)n,(r) - 2v,-(r) • S,(r) 



here and Q are auxiliary fields connected with spin and charge on the site j. 

We consider the spiral type of incommensurate magnetic order, which is the rotation of order 
parameter in the plane, modulated with some wave vector Q = {Qx,Qy) (the superposition 
states of the rotation and the ferromagnetic component perpendicular to the plane have high 
energy [14]). These states coincide with MF approximation states at zero temperature [5] 



Bj = cos (QRj) + e.y sin (QRj). 



(5) 



The two-field formalism allows us to obtain a Hartree-Fock approximation in the limit of van- 
ishing temperature [15] and phase separation [5]. From here it is more convenient for us to 
make an identical transformation and to denote 



^0 = XI ^k,tCk,t^k,t + -UMJ2 (Ck,t^k+Q,i + C] 



k+Q,i^k 



(6) 



n 



int 



Up 



/4- S, 
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with i^k,t(4.) = ^k + Unf,/2, = — 2t(cos + cos ky), where Ue is a density and M is a magneti- 
zation. Introducing Green function for Gq{z) = {z — 'Ho)~^ and making a standard calculation 
[16] we obtain 



Z*[v,C] = exp{-/3fi(v,C)}, fi(v,C) = fio - Tr In [1 - G'°^.(za;„)V,,(v, C)] 



(7) 



Here Qq = —1//3 In Tr exp [—/3(ifo ~ /^^)] is a thermodynamical potential of noninteracting 
electrons, Un = n{2n + l)//3 are fermion Matsubara frequencies, n = 0, ±1, V^^ and are 
matrices of spin variables. For simplification of the expression (4) we use the static approxima- 
tion assuming that Vj(r) and are independent on r. We neglect the thermal fluctuation of 
the charge field Q. Thus, for each configuration of the spin field Vj we set Cjiyj) equal to the 
saddle point that is given by equation 



aF(v,c) 



U [-20 



0, 



C?(v.) 



(8) 



v-Gj. We want 



We neglect longitudinal and leave only transverse fluctuations implying 
to emphasize that is the (arbitrarily chose) j-dependent unit vector. At last, we change the 
matrix expression (1— G°j V^j)~^G°j- for the average Green function G{z) (averaging is conducted 

over transverse fluctuations of auxiliary field v < ... >„ ). The self-energy Tj{z) and G{z) is found 
from the self-consistency equation treatment by a coherent-potential approximation (CPA) 



G,,{z) =< [1 - G%iz)V,,iv)]-' G%iz) >, 



E{z)). 



(9) 



At last we can write partition function in terms integration over v 

Z = j dvexp {-(3U [v^ + C\v)] } exp {-(3n{v)}, (10) 

Q{v) = Qq + ^^l^lndet[l - G{iuJn)iV{v) - - lndet[l + S(ia;„)G(iu;„)]| . 

n 

For short of the expression (10) we miss out site index j. 

Self-consistence solution of the Eqs.(8),(9),(10) allows us to calculate all magnetic properties 
of our system. 
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Fig. 1: Magnetic phase diagram for U/t = A. Bold lines denote the temperature Neel (second- 
order phase transitions), solid lines denote boundaries of phase-separated region AF-|- {Q,7r) 
(shaded area), dashed lines denote first-order phase transitions calculated without regard for 
PS. 




Results 

We have performed numerical calculations comparing thermodynamical potential of different 
magnetic phases at different n and T, solving the Eqs. (8), (9), (10) for U/t = A. Hamiltonian 
(1) have the particle-hole symmetry {n <^ 2 — n) and we can restrict ourselves to the region 
< n < 1. The magnetic phase diagram is presented in Fig. 1. At zero temperature the 



antiferromagnetic (AF) phase exist only at half-filling [n = 1) [5, 14]. The region of AF states 
is intensively grown with temperature. The dependences of on T and n are shown in the 
Fig. 2. We can see that AF — {Q,Q) phase transition is a first-order transition. From our 
calculation it is clearly that {Q,Q) — > {Q,'^) and AF — )■ (Q, vr) transitions are also first-order 
transitions. This is confirmed by the dependence of chemical potential /x on number of electrons 
n Fig. 3. We can see two regions of fi instability AF — )• {Q,Q) and {Q,Q) — {Q, vr) phases, 
which are characterized by negative derivative of /i (n = (0.92 — 0.9,0.87 — 0.86)). Therefore 
these regions must consist of two spatially separated phases. Boundaries of PS and proportion 
of phases can be obtained through the Maxwell construction. 



The Fig. 3 shows dependence of the chemical potential with account to phase separation. 



In summary, our investigations demonstrate that the magnetic phase transition qualitative 
coincides with experiments, notably the transitions from paramagnetic state to AF one and 
then to PS at low temperature [1, 2, 3, 4], see Fig. 1. However the phase separates to AF and 
{Q,7r) phases, whereas experiment shows separation to AF and {Q,Q) phases and a region 
of separation to {Q,Q) and (Q,vr) phases [1]. Behavior of the chemical potential (two regions 
of instability in Fig. 3) indicates a high possibility of the existence of PS AF +{Q,Q) and 
iQ,Q) + {Q, tt) phases at some choice of Hamiltonian parameters (1), in particular with account 
of a next-nearest-neighbor hopping [5]. 
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